pdf("Figure_8a.pdf");
res = cbind(c(log10(36142395018),log10(874707)), c(log10(53435495533),log10(1937727)), c(log10(71799469240),log10(4406657)), c(log10(90231871037),log10(8446127)), c(log10(107743851807),log10(14897016)));
colnames(res) = c("e=1","e=2","e=3", "e=4", "e=5");
par(cex=1.3);
barplot(res, main="", xlab="mrFAST", ylab = "# of Locations (Log10 Scale)", col=c("black","lightgrey"),legend = c("Seed Locations","True Locations"), ylim=c(4,15),beside=TRUE,xpd=FALSE);
dev.off();

pdf("Figure_8b.pdf");
res = cbind( c(log10(41316356102), log10(107137141), log10(874098)), 
	     c(log10(93552878043), log10(170194567), log10(1134500)),
	     c(log10(165517866791), log10(252103203), log10(4266550)),
	     c(log10(258281721852), log10(407206417), log10(8444593)),
	     c(log10(372484737520), log10(896612633), log10(14898619))

	);

colnames(res) = c("e=1","e=2","e=3", "e=4", "e=5");
par(cex=1.3);
barplot(res, main="", xlab="fastHASH with AF", ylab = "# of Locations (Log10 Scale)", col=c("black","darkgrey", "lightgrey"),legend = c("Seed Locations", "Seed Locations after AF","True Locations"), ylim=c(4,15),beside=TRUE,xpd=FALSE)
dev.off();

pdf("Figure_8c.pdf");
res = cbind( c(log10(41316356102), log10(1083533920), log10(106276000), log10(873735) ),
			 c(log10(93552878043), log10(3171096419), log10(165717191), log10(676825) ),
			 c(log10(165517866791), log10(7123623298), log10(227594084), log10(4358746) ),
			 c(log10(258281721852), log10(14131842488), log10(310217876), log10(8443682) ),
		     c(log10(372484737520), log10(25913150934), log10(437640520), log10(14897016) )

);

colnames(res) = c("e=1","e=2","e=3", "e=4", "e=5");
par(cex=1.3);
barplot(res, main="", xlab="fastHASH with AF and CKS", ylab = "# of Locations (Log10 Scale)", col=c("black",  "white", "darkgrey", "lightgrey"),legend = c("Seed Locations", "Seed Locations after CKS","Seed Locations after AF", "True Locations"), ylim=c(4,15),beside=TRUE,xpd=FALSE);
dev.off();

pdf("Figure_6.pdf");
par(cex=1.5);
res = cbind( c(4.7309027778	,3.6504065041	,4.0613496933, 6.047682119),
	     c(8.6398809524	,5.2539370079	,4.24419778, 10.712952159),
	     c(15.8816371681	,7.2114661654	,2.6049729296, 19.043441227),
	     c(11.4362969753	,2.7016089712	,1.4685575365, 13.704610951),
	     c(10.5433673469	,2.0491208485	,0.8974750782, 13.080603224)
	);

colnames(res) = c("e=1","e=2","e=3", "e=4", "e=5");
barplot(res, main="", xlab="e: edit distance", ylab = "fastHASH speedup factor over mrFAST 2.1.0.6", col=c("black","darkgrey", "lightgrey", "white"),legend = c("Dataset 1", "Dataset 2","Dataset 3", "Dataset 6"), ylim=c(0,24),beside=TRUE);
dev.off();

pdf("Figure_7.pdf", height=6, width=9);
res = cbind( c(3.8),c(3.33),c(2.32),c(3.22), c(3.12) );
colnames(res) = colnames(res) = c("mrFAST-2.5.0.0\n(fastHASH)", "mrFAST-2.1.0.6", "Bowtie-9.12.8", "BWA-0.6.1", "RazerS-3.0");
barplot(res, main="", cex.names=1.1 , cex.axis=1.1,col=c("grey"), beside=TRUE, args.legend = list(cex=1.6 , x = "topleft"), ylim=c(0,5))
title(ylab = "Memory Usage [GB]",cex.lab=1.4)
dev.off();

pdf("Figure_9.pdf");
par(cex=1.5);
res = cbind( c(1083533920 / 41316356102 * 100), 
	     c(3171096419 / 93552878043 * 100),
	     c(7123623298 / 165517866791 * 100),
	     c(14131842488 / 258281721852 * 100),
	     c(25913150934 / 372484737520 * 100)

	);

colnames(res) = c("e=1","e=2","e=3", "e=4", "e=5");
barplot(res, main="", xlab = "e: edit distance", ylab = "% of Locations", col=c("darkgrey"),legend = c("Seed Locations after CKS"), ylim=c(0,9),beside=TRUE,xpd=FALSE)

dev.off();
